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Abstract 

We have given a novel formulation of the exact solutions for the lambda, vee and cascade 
three- level systems where the Hamiltonian of each configuration is expressed in the 5(7(3) 
basis. The solutions are discussed from the perspective of the Bloch equation and the atom- 
field entanglement scenario. For the semiclassical systems, the Bloch space structure of each 
configuration is studied by solving the corresponding Bloch equation and it is shown that 
at resonance, the eight-dimensional Bloch sphere is broken up into two distinct subspaces 
due to the existence of a pair of quadratic constants. Because of the different structure of 
the Hamiltonian in the 5*17(3) basis, the non- linear constants are found to be distinct for 
different configurations. We propose a possible representation of the qutrit wave function 
and show its equivalence with the three-level system. Taking the bichromatic cavity modes 
to be in the coherent state, the amplitudes of all three quantized systems are calculated by 
developing an Euler angle based dressed state scheme. Finally following the Phoenix-Knight 
formalism, the interrelation between the atom-field entanglement and population inversion for 
all configurations is studied and the existence of collapses and revivals of two different types 
is pointed out for the equidistant cascade system in particular. 
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1. Introduction 

In the atom-field interaction scenario, the level structure of an atom leads to the prediction of 
a wide range of experimentally verifiable coherent phenomena. Probably the most notable among 
them is the observation of the collapse and revival of the Rabi oscillation [1] which unequivocally 
proves the granular structure of the photon. This phenomenon is indeed a prediction of the Jaynes- 
Cummings model - an idealized two-level system consists of an atom with two distinct quantized 
levels interacting with a monochromatic quantized cavity field [2,3]. An immediate extension of the 
two-lcvcl system is the thrcc-lcvcl system which is generally classified into lambda, vee and cascade 
configurations, respectively. Such configurations are in the purview of present studies because they 
exhibit a rich class of coherent phenomena such as, two photon coherence [4], double resonance 
process [5] , three- level super-radiance [6] , resonance Raman scattering [7] , population trapping [8] , 
tri-level echoes [9], STIRAP [10], quantum jump [11], quantum zeno effect [12], electromagnetically 
induced transparency [13] etc. Prom these studies it is quite transparently obvious that increase 
of the number of levels not only can generate a large number of quantum-optical effects, but 
also enables us to develop a suitable control mechanism which is extremely important from the 
experimental point of view. Thus the three-level configuration, the simplest representative of the 
multi- level system, demands careful inspection from time to time in its own right. 

It is well known that the Hamiltonians of the lambda, vee and cascade three-level systems 
types can be described using the atomic basis operator, a^i, = |/x >< v] (/u, v = 1,2, 3), where the 
solution is carried out with the two photon resonance and the equal detuning conditions as the 
supplementary conditions [14,15]. Apart from this treatment, another equivalent way of dealing 
with the three-level system is by the Bloch equation technique, where the eight Bloch vectors are 
defined on the eight dimensional Bloch sphere S*^ [16-18]. This method was first initiated by Eberly 
and Hioe who pointed out the relationship of the three-level system with the SU{3) group [16]. 
Their investigation revealed that the quadratic Casimir of the SU{3) group is manifested through 
the existence of some nonlinear constants, which gives rise to a non-trivial structure of the Bloch 
space of such a system [18,19]. Later, this result was obtained by solving the pseudo-spin equation 
[20] and also by the Floquet theory technique [21]. However, in the Bloch equation approach, the 
lambda, vee and cascade three-level systems are found to be generated by changing the position of 
the intermediate level E2, i.e., the energy levels are arranged as, E2 > E3 > Ei, Ei > E3 > E2 and 
E3 > E2 > El, as shown in Fig.l of Ref.[18]. In consequence, irrespectively of the configuration, 
the interaction term for any one of these thrcc-lcvcl systems in the atomic operator basis is given by, 
Hf = gi\l >< 2\+g2\2 >< 3\+h.c. {A = A,V and S). The pitfah of having identical structures of 
the Hamiltonians for different configurations is that this leads to same set of non-linear constants, 
which is undesirable because the three-level systems are intrinsically different from one another. 

Apart from that, there is another reason for studying the Bloch space structure of three-level 
systems. In quantum information theory parlance, the qubit is "designated" by various points 
on the aforesaid unit Bloch sphere [22,23]. A natural extension of the qubit is the qutrit system, 
which is generally expressed in the computational basis: ]— >, ]0 > and )-|- >; it has drawn wide 
attention in recent years [24-26]. Although there exist several suggestions for implementing the 
qutrit system either by treating it as the transverse spatial modes of single photons [27], or through 
the polarization states of the biphoton field [28,29], significant progress has been made mainly by 
identifying the qutrit with the three-level system driven by bichromatic laser fields [30] . However, 
in spite of the significant progress, a proper definition of the qutrit wavefunction and its explicit 
relation with the Bloch space structure of the three-level system are not available. 

Recently a number of quantum-optical systems have come under intense scrutiny with a view to 
the experimental implementation of various quantum information protocols where the concurrence 
is considered the most useful dynamical parameter of evolution [31]. The primary reason for 
such a study is that it is capable of deciphering the nature of the entanglement between two 
non-local objects which is the key resource in communicating information. For example, in the 
phenomenon of entanglement sudden death, the entanglement between two non-local two-level 
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systems is studied by considering the time evolution of the concurrence [32]. Since the two- 
level system is essentially an atom-field composite system, it is equally important to understand 
the dynamics of the entanglement between its constituents. As regards the two-level system, 
the correlation between the atom and the field was first considered by Gea-Banacloche in terms 
of purity [33] and then more systematically by Phoenix and Knight using the notion of entropy 
[34,35] . The latter authors also argued that the entropy can be used as an operational parameter to 
quantify the atom-field entanglement which is constrained by the Araki-Lieb limit. Since the three- 
level configuration is essentially a pump-probe system at heart, if the entanglement between two 
non-local such systems is realized experimentally, then it is possible to maneuver the entanglement 
externally. Howc!V(?r, the entanglement dynamics of two non-local three- level systems still remains 
an open issue due to the lack of availability of an appropriate definition of the concurrence of 
the three-level systems [31,36]. Therefore, prior to addressing the problem of the entanglement 
of two such non-local systems and develop a consistent theory of coherent control mechanism for 
engineering it by means some driving field, it is worth studying the entanglement scenario of the 
constituents of a single three- level system following Phoenix and Knight [34] . 

The primary objective of this paper has two parts. First, we shall study the structure of the 
Bloch space for all three-level systems while taking different Hamiltonians for different configura- 
tions. Unlike in existing treatments [16,18], it is explicitly shown that if the Hamiltonians of the 
three-level systems are expressed in the SU{3) basis with the same energy condition, this leads 
to distinct non-linear constants for different systems. We also give a convenient representation 
of the qutrit system and discuss its properties. Secondly, after developing an Eulcr matrix based 
dressed state formalism, we discuss a systematic scheme for calculating the atomic entropy follow- 
ing Phoenix-Knight [34] and show how the increase of the number of levels influences the atom-field 
entanglement and the population inversion for all three configurations. 

The remainder of the paper is organized as follows. In Section 2 we reconsider the Hamilto- 
nians of the semiclassical lambda, vee and cascade three-level system types and their quantized 
versions expressed in the SU{3) basis. Section 3 obtains the solution of the Bloch equations for the 
semiclassical configurations in order to provide various non-linear constants and gives a possible 
representation of the wave function of the qutrit system. In Section 4, we proceed to obtain the 
quantized configurations by introducing an Euler angle based dressed state formalism and compare 
the corresponding atomic entropy with the population inversion for various initial conditions. In 
Section 5, we numerically study the atomic entropy with the population inversion for various initial 
conditions and give a possible estimate of the collapse and revival times for the cascade system in 
the high field approximation. Finally we summarize the main results of the paper and discuss the 
outlook. 

2. The models 

The time-dependent Hamiltonian H{t) can be expanded in the basis of the Lie operators as 

N 

H(t) = ^ hi(t)Li, where hi(t) are the linearly independent complex valued functions and the Li 

i=l 

are the generators which satisfy the Lie algebra [Li, Lj] = iwijkLk- The algebraic structure of the 
Hamiltonian ensures the existence of a unitary operator U{t) = YiiLi Gxp['i9i{t)Li], where gi{t) is 
the scalar function which is to be evaluated for any specific model. Thus the time-dependent wave 
function can be evaluated either by the unitary method, i/j{t) = U{t)ip{0) [37] or by using suitable 
dressed state formalism [38,39]. In this Section, we develop the Hamiltonian of the semiclassical 
and quantized three-level systems where the atomic levels are expressed in the SU{3) basis. 
It is customary to consider the three-level system as two two-level systems in the dipole approxi- 
mation and the Hamiltonian of the semiclassical lambda system is given by 

H^ = hiiOiVa + W2T3) + di3 • Ei(Oi) + das • £3(02), (1) 
where, dis = (V+-|-V_) and d23 = (T+-|-T_) be the dipole operators representing the transitions 
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1 O 3 and 2 o 3, and Ei(r2i) — Ei+ + Ei_ and ^2(^2) = E2+ + E2_ are the bichromatic classical 
electric fields with mode frequencies ili and 172, respectively. Similarly we have 



= hicoiVa + W2C/3) + di3 • Ei(ni) + di2 • E2(r!2), 



(2) 



for the vee system, where dis = (V_|-+V_), di2 — (U++U„) be the dipole operators representing 
transitions 1 4-)- 3 and 1 O 2, respectively. Finally for the cascade system we have 



H= = n(cJi J73 + W2T3) + di2 • El ) + d23 • E2 (^2 



(3) 



where di2 — (U+ + U_) and d23 = (T+ + T_) are 1 o 2 and 2 3 transitions, respectively. 
In Eqs.(l-3), huji{^ —E^),hw2(= —E2) and h{(jj2 + wi)(= -E3 ) are the energies of the three 
levels of the lambda system, 1twi{= EY), fiu;2{= E2) and —h{u;2 + uJi){= E^) are those of the 
vee system and lluji{= —Ef),h{ui — W2)(= Ef) and hi02{= Ef) are the energies of the cascade 
system, shown in Fig.l. Here we note that the energy levels of all three configurations satisfy 
an unique condition Ei < E2 < E3, and this leads to three distinct interaction Hamiltonians, 
namely, = gi\l >< 3| + g2|3 >< 2| + /i.e., iJf = gi\l >< 2| + 52II >< 3| + h.c. and 
Hf = gi\l >< 2\+ g2\2 >< 3| + h.c, respectively. These interaction Hamiltonians in the atomic 
operator basis are precisely identical to the dipole interaction terms in Eqs.(l-3) which are expressed 
in the SU{3) basis. In Eqs.(l-3), we have defined the SU{3) shift vectors [40], 



T± = i(Ai±iA2), 



f/± = i(A6±iA7), 



T3 = A3, 



V3 = (^/3A8 + A3)/2 



F± = i(A4±iA5), 

Us = {VSXs - A3)/2, 



where the Gell-Mann matrices are given by. 



Ai 



A4 = 



1 

10 



1 



1 



A5 = 



-i 

1 


-/■ 



1 



A3 



A6 = 



1 
0-10 





1 

10 



(4) 



A7 = 




-i 
i 



1 

V3 



1 
10 
0-2 



These matrices follow the commutation and the anti-commutation relations 



[Ai, Aj] — 2i/jjfeAfe, 



{Ai, Xj} — + 2dijkXk, 



(5) 



(6) 



where dijk and fijk = 1, 2, .., 8) represent completely symmetric and completely antisymmetric 

structure constants, respectively, which characterize SU (3) group. 

In the rotating wave approximation (RWA) , the Hamiltonians of semiclassical system can be 
written as [39], 



= hini -uji- L02)Va + h{n2 - wi - w2)T3 + HA^Vs + a^Ts) 



+hKiV+ exp(— iQit) + hK2T+ exp{—ifl2t) + h.c. 



for the lambda system, 



1 - wi - W2)V3 + h{^2 - - a;2)U3 + h{AXVs + A^C/3) 
+?iKiV+ exp(— iOif) + hK2U+ exp(— if22t) + h.c. , 



(7) 



(8) 
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for the vee system, and 

= - wi + LJ2)U3 + h{Q.2 + wi - a;2)T3 + ;i(Af J/g + AfTg) 

+7lKif7+ exp(— zfiif) + fiK2T+ ex.p{—ifl2t) + h.c. , (9) 

for the cascade system, respectively. In Eq.(7-9), Kp (p = 1,2) are the couphng parameters and 
Af = (2wi +UI2 - ^1), A^ = (wi + 2w2 - ^2) for a = A,V and Af = (2wi - a;2 - l^i) and 
Af = (2w2 — wi — $^2) are the respective detuning frequencies. 

For the quantized system, the Hamiltonian of the lambda system is given by 

H^ = H^ + (10) 

where, 

2 

Hq = h{ni — coi — C02)V3 + h{fl2 — wi — <jJ2)T3 + ^ ^a^Oj (llo) 

Hf- = n{A^V3 + A^Ta) + h{giV+ai + 52^+02) + h.c, {lib) 

with aj and aj {j = 1,2) the creation and annihilation operators of two cavity modes, gj are the 

coupling constants and D,j arc the mode frequencies. 

Similarly the Hamiltonian of the vee system is given by 

= + hY (12) 

where, 

2 

Hq = h{Qi - coi- W2)V3 + h{Q2 - wi - ^2)^/3 + ^2 ^j^l'^j (13a) 

i=i 

hY = HAYV3 + A^Us) + h{giV+ai + g2U+a2) + h.c. (136) 
Finally the Hamiltonian of cascade system in the same approximation is given by 

H^ = Hf + Hf (14) 

where, 

2 

Hf = h{fli +UJ2- i^i)U3 + h{Q.2 +0J1- 0J2)T3 + ^ ^ja]aj (15a) 

Hf = h{AfU3 + AfTs) + h{giU+ai + 52^+02) + h.c. (156) 

Using the algebra of Gell-Mann matrices it is easy to see that for the lambda and vee systems 
(a = A and V) 

[HS,Hf]=0, (16) 
provided the two-photon resonance condition, namely, 

A? = -Ag, (17) 
is satisfied. Similarly for the cascade system we must have 

[Hi,Hf]=0, (18) 
which requires the equal detuning condition, i.e.. 
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(19) 



Thus the two-photon condition and the equal detuning condition ensure the commutativity of 
the free and interaction parts of the Hamiltonian showing that the quantized models are exactly 
solvable. We now proceed to discuss the Bloch space structiire of all thrcx^-k^vcil configurations 
by exploring the possible non-linear constants obtained from the semiclassical Haniiltonians in 
Eqs.(7-9) written in the SU{3) basis with a unique energy condition, Ei < E2 < E3. 

3. Bloch equation and non-linear constants 

Let the solution of the Schrodinger equation of the semiclassical three-level system described 
by the Hamiltonians in Eq.(7-9) be given by, 



^^{t) ^ C^(t) I-) + C^{t) |0) + C^{t) 1+) 



where C^,Cq and {A = A,V and S) are the amplitudes of the atomic states which satisfy 
the normalization condition |C^p -|- \Cq\'^ + |C^P = 1. In Eq.(20) the basis states of the lower, 
middle and upper atomic states are given by. 



>= 




|0>= 




+ >= 




(21) 



The exact evaluation of the probability amplitudes enables us to calculate the density matrix 

defined as 

p^{t) = |*^(i) >®< ^^{t)\. (22) 

To start with, we consider the dressed wave function of the lambda system obtained by a 
unitary transformation 

^^{t) = uim^iO), (23) 
where the unitary operator is given by 



UA{t) = exp[--{{2n2 - ni)T3 + (20i - n2)V3)t]. 
The time-independent Hamiltonian of the lambda system described by Eq.(7) is found to be 

ifA(o) = -huiiiA + uiH^{t)UA 



(24) 



h{A^ + A^) hK2 



hK2 
flKl 











(25) 



Similarly, the unitary operator of the vee system described by the Hamiltonian in Eq.(8) is, 

Uv{t) = exphf ((2r!2 - ^i)U3 + (20i - 02)V^3)t], (26) 

and the corresponding time-independent Hamiltonian is given by 

^h{2AY - AX) hKi 

|?i(2A^ - AY) hK2 

hKi hK2 -\h{AX + AX) 



HV{0) = 



(27) 
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Also the Hamiltonian of the cascade system in Eq.(9) can be made time- independent by using the 
transformation operator 

Usit) = exp[~{{ni + 202)73 + (2fii + fl2)U3)t\. (28) 



and the corresponding Hamiltonian is 

" U(Af + 2Af) 



iJ^(O) = 



hK2 





hK2 

h{Af - Af ) 







hKi -lh{2Af + Af) 



(29) 



Thus we have three distinct Hamiltonians for three different configurations. 
To obtain the Bloch equation, we define the generic SU (3) Bloch vectors 

Sf{t) = Tr[p^m], (30) 

where p"^ be the density matrix given by Eq.(22) which satisfies the Liouville equation, 



Eq. (30) can be equivalently expressed in terms of a Bloch vector. 



P^(i) = ^(l + |ES^(W- 



(31) 



(32) 



i=l 



Substituting Eq.(32) in Eq.(31) and making use of the Hamiltonian of the requisite model we 
obtain the Bloch equation 



dt 



(33) 



where M^^ be the eight-dimensional anti-asymmetric matrix. For the Hamiltonian of the lambda 
system given by Eq.(25), the matrix reads 






A'^ 
















-Kl 





-A^ 





2k2 








-Kl 













-2k,2 

































AA 







K2 











Kl 


-AA 





-K2 







■\/3ki 





Kl 








K2 







(A^-AA) 





Kl 








-K2 





-(AA- 


A^) 




















— V^Ki 














(34) 



Similarly, from the Hamiltonian in Eq.(27) we find the Bloch matrix 





(Ar - Ar) 








-K2 





-Kl 





(Ar - Ar) 








K2 





-Kl 




















-Kl 





K2 
















Ar 











K2 





Kl 


-Ar 
















Kl 














A^ 
^2 





Kl 





-K2 








^2 





\/3k2 














-\/3ki 





-\/3K2 






(35) 
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for the vee system, and from Eq.(29) we get 






Af 








-Ki 











-Af 
























-2k2 














Kl 
















(Af + Af ) 





K2 





Ki 








-(Af +Af) 

























K2 





Af 













-«2 





-Af 

























— \/3ki 






(36) 



for the cascade system. 

The algebraic structure of SU{3) group allows the existence of a set of quadratic Casimirs which 
will appear in the form of quadratic constants. To construct them from eight Bloch vectors, we 
have to search through a total of (g_^)!„! combinations which will form a tuple. However, because 
of the large number of such combinations, finding the exact number of such tuple by solving the 
Bloch equations (34-36) is quite onerous. Therefore, to determine all possible non-linear quadratic 
constants, we have developed a Mathematica program to carry out an extensive search and obtain 
the following results for n = 3 and n = 5 only. 

At resonance (A^|^ = = A2 ), for the lambda system the Bloch space is constituted of two 
parts, one being of two-sphere 5^, 

Sf{t) + Sf{t) + Sf{t) = Sf{0) + Sf{0) + Sf{0), (37a) 
and the other the four-sphere 5^, 

sfit) + sfit) + Sf{t) + Sf{t) + S^\t) = 

S^\0) + Sf{0) + St\0) + S^\0) + S^\0), (376) 

respectively, where S^-{0) are the Bloch vectors at t = which are to be evaluated in terms 
of the probability amplitudes. On noting the fact that the density matrix can be written as 
p^{t) = W{t)p^{0)U{t), the Bloch vector in Eq.(30) becomes, 

5f (0) = rr[p^(0)A,]. (38) 

Inserting Eq.(5) and the density matrix from Eq.(22) into Eq.(38) at t = 0, all constants Sf"{0) in 
Eq.(37) can be expressed in terms of probability amplitudes, 

Sf + Sf + 5^' = iC^{0fC^{0f + 4(7^(0)2(7^(0)2; (39a) 

and 

^2''" + + + S^' + S^' = ^(C^(0)2 + C^{0f + c^ioff 

-3C^{0fc^{0f - 3C^{0fcl{0f, (396) 
respectively. Similarly at resonance (Aj^ = = A2'), for the vee system we have two Bloch spheres, 

SY' + Sf + = ^C^iOfcXiOf + 4CX{0fcY{0f (40a) 

and 

sf + sf + sf + sf + sf = ^(cr (0)2 + C^iOf + CliQff 

-3C^{0fc+^{0f - 3(7^(0)^(7^(0)'. (406) 
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Also at resonance (Af = = Af), for the cascade system we have 

Sf + Sf + Sf = 4Cm'Cf{0f + 4Cf{0fCl{0f (41a) 

and 

sf + sf + sf + sf + sf = ^{C^iOf + Cf (0)2 + c^ioff 

-SCm'ciiof - 3Ci{0fcf{0)\ (Alb) 

From Eqs. (39-41) it is evident that at resonance (A^ = = A^), the seven-sphere 5'' is broken into 
two-subspaces, namely, two-sphere .S^ and four-sphere <S*, repectively, each with distinct norm. 
We also note that for 3- and 5-tuple of the Bloch vectors constituting such two quadratic constants 
are in contrast to the case for [16-20], different for different three- level systems. Finally, we note 
that for the off-resonance case (Af ^ ^ ^2)^ solutions of the Bloch equation of all three 
configurations satisfy, 

sf + sf + sf + sf + sf + sf + sf + sf = ^{C^iOf + C^iOf + (7^(0)2)2. (42) 

The normalization condition, {0)^ + Cq {0)"^ + C^{0)'^ = 1, readily gives that the radius-squared 
of the eight-dimensional Bloch sphere is | , regardless of the configuration. 

3.1. Qutrit wave function 

Having knowledge of the Bloch space of the three-level system, we now consider the qutrit wave 
function which is the qubit counterpart of the two-level system, we define, in the three dimensional 
Hilbert space H^, the four-parameter normalized qutrit wavefunction, 

>= cos ||- > + sin I sin | sin |e^*|0 > + 

(sm y c;os Y + « sm — sm — cos y )|-|- >, (43) 



where ^1 and 62 lie between to tt with an arbitrary value of the phase angle (f>. From Eq.(43) 
it is evident that, the pure qutrit states |— >, |0 > and |-|- > can be obtained by taking the 
coordinates to be (0,^1,^2), (tt, tt, tt) and (77,0,^2), respectively. This indicates that for a qutrit 
system, the state |— > may be any point on the (^1,^2) plane, |0 > corresponds to a well-defined 
point and |-|- > lies on a line for any value of ^2- Furthermore, it is worth noting that by taking the 
angles (^1, ^2) to be (tt, tt), the qutrit wave function in Eq.(43) is reduced to a qubit wavefunction, 
which indicates the consistency of the qutrit wave function. Using Eq.(22), various elements of the 
density matrix of the qutrit wave function are given by 

11 2 ^1 
Pt =cos —, 

22 . 2 ^0 . 2 ^1 • 2 ^2 
prp = sm — sm — sm — , 
'^^ 2 2 2 

1 6 

- (3 -|- cos ^1 -|- cos 62 — cos 9i cos 62) sin^ — , 

12 21* ^ iHt ■ n . ^1 . ^2 

Pj. = Pj. = -e'*" sm 9q sm y sm y , 

-id> ■ 2 . Oi Oi . 9i O2. . 02 
e sm -— sm -— (cos -z- +ism — cos -— ) sm -— , 



p¥ = 
pf = pf* = 
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Pj, = Pj, = - sm^o(cos y + ^ sin y cos y). (44) 

It is worth noting that this wavefunction is normaUzed, i.e., Tr[pT] = 1, which also satisfies the 
pure state condition p^ = p^. Finally plucking pr back in Eq.(30) we obtain, 

Sf + Sf + Sf + Sf + Sf + Sf + Sf + Sf = I (45) 

which is precisely same as Eq.(42) obtained by solving Bloch equation for different configurations 
of thrcc-lcvcl systems. This establishes the complete equivalence between the qutrit with the 

thrcc-lcvcl system. 

4. The Euler angle approach to dressed state formalism 

At t = 0, the atom- field wave function of the quantized three-level system is factorized as 

oo 

Vaf(0) = (C^(0)|->+Co^(0)|0>+C^(0)|+>)0 CnC^\n,m>, (46) 

n,m=0 

where Cf be the atomic amplitude with level indices i = — , 0, + and Cm,Cn are the amplitudes 
of the bichromatic fields which are in the coherent state with field modes m and n, respectively. 
With time, the atom and field will be entangled with each other and the instantaneous atom-field 
wavefunction for i > is given by 

oo 

Mt)= Yl {<:?™:^„WKn,->+C^„(t)|m,n,0>+C+;^„(t)|m,n,+ >}, (47) 

n,m=0 

where C^^„(t) are the amplitude in the interaction picture which are to be evaluated. Following 
the Phoenix- Knight formalism [34] , the atomic entropy of the three- level system is given by 

SA{t) = -Tr[pA{t)lnpA{t)], (48) 

where pA be the reduced density matrix of the atomic system obtained by taking trace over the 
two cavity modes: 

PA{t)=TrF„FdpAF{t)]. (49) 

In Eq.(49) the density matrix of the atom-field system is given by, 

PAFit) = \^AF{t) >(S>< iJAFit)\. (50) 

Finally the atomic population inversion between 1 -H- 2, 2 -o- 3 and 1 3 levels are given by 

=< Ts >= Trp^F, [pAF{t)T3], (51a) 

W^i =< Us >= Trp^F, [pAFms], (516) 

W^s =< Vs >= Trp^F, [pAFim]. (51c) 

In Appendix we have developed an Euler matrix based dressed state formalism to calculate the am- 
plitudes for three possible choices of the basis states shown in Fig.la-c. In the following subsections 
we present their explicit derivations for all three configurations at resonance (A^ = = A^). 
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4.1. The Lambda system 



We first consider the interaction Hamiltonian of the lambda system given by Eq.(ll) for all 
three possible initial conditions discussed in the Appendix. For the initial condition given by 
Eq.(A.la), the Hamiltonian in the number state representation at resonance {A'l = = A2) is 
given by 

g2y/n+l giy/m 

Ht'^ = 52\/nTT . (52) 

gwpm 



The eigenvalues of the Hamiltonian are A^^ = —\Jg\m + g\{n + 1)(= — O^), A^2 = ^'^^ A^ 
g\m + g\(n + 1)(= f2^) and the various elements of the orthogonal matrix in Eq.(A.8) are 

A 



a 



11 



1 

0, 



_ A 
12 



92 



K+1 

2 ■ 



*31 



1 

V2' 



"22 =Wt 



n+l 
2 ' 



_ A 
"13 = 


31 


fin 
1 2 ' 


A 

"23 = 


92 




_ A 
"33 = 


91 


pfn 


1W\ 


2 • 



(53) 



The corresponding Euler angles are found to be 



51 



= arccos[^^/ — ], = - arccos[-52^ 



n+l 



0^'+5i(n+l) 



^^3 = arccos[— 



2n + 2 



+9i(n + l) 



(54) 



Thus, when the system is initially in the lower level, the corresponding time-dependent amplitudes 
in Eq.(A.2a) are given by 

,_A gi(n + l) + ff>cos»^^ 

am,n(*) = ' (55a) 



'^m-l.n+lW 



g\gi\Jfn{n + l)(cosrif^f — 1) 



'^m-l.nW 



.gi^rn 



(556) 
(55c) 



For the initial condition given by Eq.(A.lb), the interaction Hamiltonian in the number state 
basis is 

52 giy^m+ 1 

g2V^ . (56) 



giVtn + 1 







The eigenvalues of this Hamiltonian are given by A^ 



01 



Vdii^ + 1) + 



-r!^), =Oand 



Aq3 = ^y gf{m + 1) + g2n{= fig ) the various elements of the orthogonal matrix of Eq.(A.8) 
are given by 



a 



11 



^21 



*31 



1 
0, 



'12 



32 fn 
Tt^V 2' 



1 

V2' 



a 



32 



92 fn 
V 2' 



*13 



91 



m+1 
2 ' 



*23 



(57) 



"33 = 



_ 91 



1 / "i+i 
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which readily give the corresponding Euler's angles; 



91 m+l 



'oi = arccos[— ^W— — ], = - arccos[-£?2. 



], 6*^3 = arccos[-5i2. 



2n 



(58) 

Thus, when the system is initially in the middle level, the corresponding time-dependent amplitudes 
in Eq.(A.2b) are given by 



'^rn+l.n-lW — 



<n(i) = 



9i92^{m + l)n(cos n^t - 1) 



gf{m + 1) + cos fig t 



.g2\/nsmSlQt 



OA 



Finally the interaction Hamiltonian for the initial condition Eq.(A.lc) is given by 





g2Vn + 1 gWm + l 
g2 VnTT 
gWm + 1 



(59a) 

(596) 
(59c) 

(60) 



and the eigenvalues of the Hamiltonian arc given by = —\/g\(jn + 1) + g^in + !)(= — O^), 
\%2 = and X{l_3 = ^ygl{m + 1) + gi{n + 1){= 0|) with various elements of the orthogonal matrix 
are 



'11 



1 

V2' 



a. 



ati^ = 0, 
+ A 

"31 =-- 



_ 92 / n+1 
'12 - Y 2 ' 

0^22 = TP^VmTT, 



*13 



^23 



91 / m+l 



a 



A 



32 - n5 



(61) 



92 /n+1 



Then the corresponding Euler angles are found to be 



^33 - m 



gi I m+l 



^+1 



arccos 



0| 



51 m + l A r 

-z—], 6'+2 = -arccos[-5i2. 



n + 1 



?+3 = arccos [-32, 



2n + 2 



of +5|(n + l)^ 



(62) 



Thus, when the system is initially in the upper level, the time-dependent amplitudes in Eq.(A.2c) 
are given by 



,+A . . _ _ . ffi\/m+ IsinO^t 

"•m+l,n\'') ~ * qA ' 



'^m,n+l(^) ^ 



,(j2\/>i + IsiuOi^t 



0| 



d+A„(t)=cos0^^i. 
Substituting Eqs.(55), (59) and (63) in Eq.(A.2) we obtain. 



(63a) 

(636) 
(63c) 
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oo 

E 

m=0,n=0 



+(7^ cos n^t}|m, n, + > -iC^{^i^^ ^ + 1, n, + >+^^^^ ^ ^|m, n+1, >}} 



y. CM '''' ^ 'l^'" " ^ V - 1, n + 1, >-,gL^^^" "-^ 



m=l,n=0 



QA 



m— 1, n, + >} 



, V" n n r 9192 \/n{m + l) {cos n^t- 1) .^zV^sinf^^f 
+C^' 2^ C„C„{ p |to + l,n- 1,- >-i —j^ — ^|m,n-l,+ >}. 

m=0.ri=l "O 

(64) 

Finally rearranging the summation, the atom-field wavefunction for the lambda system in the 
entangled basis is given by 

oo 

^AF{t)= E {O*)K"^'->+CrnWK^'0>+C'^™(i)Km,+ >}, (65) 

n,m— 

where the time-dependent amplitudes are given by, 

r-^ m _ r^r r 9l{n + I) + glm cos nH ^ 

W,mW — ^-^m^n ^^^^"2 T 



A gifl2Vrn(n+ l)(cos»^t- 1) ^ 

'-^m-lL'n+l ^^^^"2 ^ 



gxy/msm9.H 

-iC'^Cm-iCn ^ , (66a) 



r^oA _ ^ giff2Vn(rn+l)(cosl^^t-l) 

^n,m\'^J — <-^-'-^TO+l<-^n-l ^i^T-g 



fff(m -H) + g|n cos 

T<-^0 '^rn'^n Tt 

QA 

(666) 



_LA . X A ai \/w + 1 sin fiAt 



„ g2\/n + lsinn^t 
+C|'^C„C„cosO+i. (66c) 
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4.2. The Vee system 



For the vee system at resonance (A]^ = = ^2^), the interaction Hamiltonian in Eq.(12) with 
the initial condition Eq.(A.la), in the number state representation is given by 



H 



-V 



giVrn 

52V^ 

gi\/m gisfn 



(67). 



The eigenvalues of the Hamiltonian are given by Al^^ = —\/g\m + g\n(= — ^il^), Al^2 = and 
A1^3 = ^J'g\m^^g^(= and the mixing angles are given by 



_i = --, e_2 = axccos[--^y/n\, 6'_3 = --. 



(68) 



Thus, when the vee system is initially in the lower level, the time-dependent amplitudes are given 
by 



.gi\/rns'm ilXt 



OX 



For the initial condition Eq.(A.lb), the interaction Hamiltonian is given by 



(69a) 
(696) 

(69c) 



gi ym 

g2Vn + 1 



giVrfi gi\/n-\-\ 



(70) 



The eigenvalues of the Hamiltonian are given by Aji = —\Jg\fn + g\(n + 1)(= A^ = and 

Aq^ = \J g\m + g\{n + 1)(= Og^) and the corresponding Euler angles are 



C = - ^ . ^02 = arccos[- ^ Vn + l] , Q^z 



2' 



(71) 



Thus, when the vee system is initially in the middle level, the time-dependent amplitudes are given 

by 

MV u\ _ ,g2\/n+lsmQ,Xt 



'm,n+l 



it) 



gfm + g2{n + l) cos ClQt 



,0V _ gig2Vm{n + l){-l + cosnXt) 

"m-l,n+ll'^J — 



(72a) 
(72&) 

(72c) 



Finally the interaction Hamiltonian in the number state representation with the initial condition 
Eq.(A.lc) is given by 

gWm + 1 

H+" = g2V^ (73). 



gi^m + 1 g2^/n 
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The eigenvalues of the Hamiltonian are given by = —\J(j\{m, + 1) + g\n{= — ri!5j^), A^2 = 
and = sj g\(rn^V)^ g\n = ip^^) and the Euler angles are given by 

Q\^=-\, ^]^2 = arccos[-||^A^, ^\z = -\ (74) 

and the corresponding time-dependent amplitudes when the system is initially in the upper level 
are given by 

<';.,„(') (75.) 



A+v (+\ _ gmVm + ly/n{-l + COS Q^t) 



(756) 



Ut) = ^^^^^2 -■ (75c) 

Plugging Eqs.(69), (72) and (75) back in Eq.(A.12), the atom-field wavefunction in the entangled 
basis is given by 

oo 

^AF{t)= E {C-);,(t)|n,m,- >+C^V„ (01^,^,0 >+C+);(i)|n,m,+ >}, (76) 
where Ci^"''"^{t) be the required time-dependent normalized constants which are given by 

Crn,nV) = C-CmCn COsQ^t - lCoC„C„_l — p 

giV^sinOrt 

-ZG+Gm_iC„ —y , (''«) 



+GoOtoG„ — 

"0 



, ^ ^ ^ gig2V'm(n-H)(cosn^t- 1) 

"0 

r^-FV _ -r^ n 9Wrn + Isinnlt flig2\/^\/m + l(cos»^t - 1) 

'~'m,nV') — ~*G_Oto-|-1G„ —^7 h OoOto-|-1G„_i 

(77c) 

4.2. The Cascade system 

For the cascade system, particularly for its equidistant configuration, we have only single quan- 
tum number with m = n and gi = g2 = g. At resonance (Af = = Af ) with the initial condition 
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Eq.(A.la), the interaction Hamiltonian in Eq.(14) of the equidistant cascade system in the number 
state representation is given by 



gVn^ 

gVn- 1 

gy/n 



(78) 



The eigenvalues of the Hamiltonian are given by Ar:i = —gV2n — 1(= —fir:), X-2 = ^^id Ar:3 = 
g^/2n — 1(= i^-) and the corresponding Euler angles, which make all states entangled, 

by 



-3 

are given 



n 



9~i = — arccosfw 1, 0^^ = — arccos[— 

O'z-i = — arccosf— 



2n- 1 
3n-V 



l2n — 2 

3 = - arccos[- ^— ^] . (79) 

Thus, when the cascade system is initially in the lower level, the time-dependent amplitudes are 
given by. 



q^in — I) + q^n cos Q-~t 

dn~{t) = — '-^2 

,_H u\ .gV^sinn^t 
dn-i{t) = -r ^ , 

,-H i.\ - gV(» -!)»(-! + COS ngt) 



(80a) 
(806) 
(80c) 



Hf = 



Similarly for the initial condition Eq.(A.lb), the interaction Hamiltonian is given by 

g^ " 

gy/n gVnTl . (81) 

gVnTT J 

The corresponding eigenvalues of this Hamiltonian are Ao_ = —g^/2r^+T{= — f^o )' ^02 = and 
= g\J2n + 1(= Og ) the Euler angles are found to be 



n + 1 



'01 = -arccos[^^^], = - arccos[-^^^], 



2n+ 1, 



6*03 = — arccos[— 

Thus, when the cascade system is initially i 
given by 



'^]. (82) 
' in the middle level, the time-dependent amplitudes are 

,n- / ^ + 1 sin riS't 

(^n+iW = -i— (83a) 

"0 

d°=(f) = cosfift, (836) 

er(i) = -^^^=^. (83c) 
"0 

Finally for the initial condition Eq.(A.lc), the Hamiltonian in the number 
is found to be 



sVnTT 
5\/n+ 1 gVn + 2 
s'\/n + 2 



state representation 
(84)) 
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The eigenvalues of the Hamiltonian are given by X^i = —g\/2n + 3(= — A^2 = ^^'^ ^+3 
gy/2n + 3(= fi^) and the corresponding mixing angles are 



C,i = — arccosU/ , , o = — arccos — / 



2n + 2 

= - arccos[- ^ ^-^] , (85) 
and the time-dependent amplitudes when the system is initially in the upper level is given by, 



.+B _ ffV(n+l)(n + 2)(cos»ft-l) 



(866) 



,+H.,x g"(n + 2) +ff^(n+ 1) COS f^l^ 

< (i) = . (86c) 

Plugginging Eqs.(80), (83) and (86) back in Eq.(A.14) we obtain the atom-field wavefunction for 
the cascade system in the entangled basis; 

oo 

^Fit) = ^{C-=(t)|n,- > +Cf (0|n,0 > +C+^{t)\n,+ >}, (87) 

n=0 

where the time-dependent normalized constants are given by, 

W — *-^-<^n ^^^2 



,^E^E gVn(n-l)(cosOgt-l) 



-^C+C„_r^^^^^, (886) 

"0 



^^.^Hg> + 2)+gV+l)cosOft_ ^^^^^ 
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5. Numerical results 



To understand the relationship between the atom-field entanglement and the population inver- 
sion of the quantized lambda, vee and cascade systems, we consider the following initial conditions 
which are relevant from the experimental point of view: a) Case I: C^(0) = 1, Cq{0) = 0, C^{0) = 

0, i.e., when the system is initially in the lower level; b) Case II: (0) = 0, Cq{0) = 1, Cf{0) = 0, 

1. e., when the system is initially in the middle level; c) Case III: (0) = 0, Cq{0) = 0, Cf{0) = 1, 
i.e., when the system is initially in the upper level. The cavity modes are taken to be in the 
coherent states which are characterized by the Poisson distribution, namely. Cm = e"'"""' ^^^^r 

and Cn ~ e~'""'^^^^%ri where jofmP and |q!„P be the average photon number of two modes of 
the bichromatic fields. Finally, inserting the wavefunctions given by Eqs.(65), (76) and (87) in 
Eqs. (49-50), and then taking trace over the bi-chromatic field modes, we obtain the atomic en- 
tropy from Eq.(48) and population inversion from Eq.(51). For all configurations, the entropy lies 
between < §a ^ log2'^ and below we consider their time evolution and compare them with the 
corresponding population inversion for different cases. 

(a) The Lambda and vee system: Due to the existence of two modes in the Rabi frequency 
of the lambda and vee systems, the calculation of the collapse and revival times is non-trivial. 
Therefore, for these systems we give a numerical solution of the atomic entropy and population 
inversion only. Figs. 2 and 3 give the plots of the atomic entropies and popiilation inversions of the 
lambda and vee systems respectively for Cases I-III. From these graphs it is evident that during the 
period of collapse, the entropy reaches a minimum value which is sufficiently away from zero. This 
indicates that the system does not return to the pure state during the quiescent period of collapse. 
Furthermore, it is noteworthy from the population inversion plots that, during collapse period, the 
difference in population of the involved levels is no longer zero, but maintains a constant value. 
Both of these features are in sharp contrast to the two-level system, where the atomic entanglement 
goes to almost zero during the period of collapse and two levels are equally populated, which leads 
to a vanishing value of the population inversion [34]. 

(b) The Cascade system: Figs. 4-6 show the numerical plots of the entropy of the equidis- 
tant cascade system for aforementioned cases and compare them with corresponding population 
inversions. We note that the plot of Case-I in Fig. 4 is almost identical to that for Case-Ill in 
Fig. 6 showing that there exists an inversion symmetry for this system. In addition. Fig. 4a, b and 
Fig. 6a, b show that the system undergoes collapses and revivals of two distinct types with different 
half-width; this is not exhibited by the two level system. We now analyze the derivation of the 
enveloping function of the cascade system. 

The presence of single mode in the Rabi frequency of the cascade system enables us to get 
a reasonable estimate of the collapse and revival time in the high field approximation. In this 
approximation the mode frequency n is expanded around the average photon number n which is 
assumed to be very large, i.e., n{= |ap) >> 1 [41]. Below, we have given their detail derivation 
for Case-I and we quote the results for other cases. 

Case I: The wave function of the cascade system for this case can be easily read off from Eq.(87); 

VIfW = V'-(i)|- > (i)|0 > +V?(i)l+ >, (89) 
where the field components obtained by rearranging the cavity mode are, 

^_{t) = 2_^C„ —2 \n>, (90o) 

n=0 

m) = -^Y.^n ~ \n-l>. (906) 

x—^ q^Jnin — l)(cos fl~t — 1) , 
^~it)^Y.^n^^ -^2 = -\n-2>. (90c) 

n=0 
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Substituting Eq.(89) in Eq.(51a), the population inversion of the cascade system is, 

CO 

Wr2it) = Yl '^niPn + QnCOSn- + RnC0s2n~ ), (91) 



n=0 



where, P„ = ^ 2(2k-i7^^'' > ^ "2n-^)^^^ ^ ~ 2(2w-i)^ ' respectively. To obtain the summation 
over the cavity modes in the aforesaid high field approximation, the Rabi frequencies in Eq.(91) 
can be expanded as, 

nlc.'±^l±^ + 0{{n-nn (92a) 
\/2n — 1 

n=.^i!L+^+C)((n-^)^), (92/,) 
V2n + 1 

^%^'-^^^^+0{{n-nf). (92c) 

Vo + 2n 

In the limit h » 1, the polynomials in Eq.(91) become constant, i.e., 

113 

Finally, using the value of C„ and inserting Eqs.(92) and (93) in Eq.(91), the summation over the 
cavity mode gives the population inversion between various levels for Case-I, 

^^(->^ ^ 8 2 ^ V2n - 1 V2ri- 1 



3 _n{i-cos[-|«i_]} ^ r2fl'(n-l)i , _ . r 2s(f 
e iv^J^ X cos[^==^ +nsm[-===]], 

o v2n — 1 V 2n — 1 

W^2lr-^ W = 4 + ie-^{^-°«[7feI} X cos[^i^S=i}^ + nsin[ 



(94a) 



gt 



3 _}j{i_cos[^|ii_]} r2c/(n-l)f , _ . r 2fifi 

e xcos[ ^ — - +nsm[ — J ], 

o \/2n — 1 v2n — 1 



(946) 



W^.V)W = e-^-™^'v*=Tl} , cos[^|^ +^sin[_iL=]]. (94c) 
^ ' yzn — 1 ^/2n — 1 

Prom the two enveloping functions given by Eqs.(94a) and (94b), we note that there are two distinct 
times of collapse and times of revival. In the limit n >> 1, the two times of collapse are 



V2(2n- 1) ^ 2 
gVn fit' 



!liz) = ^ 1, (956) 



and similarly, two times of revival are 



2Vm^ 2V2fi 

t-,(\ = TT ~ TT, (96a) 

^'99 

H-) = ^ = ^ - — -TT. (966) 
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In contrast, the population inversion between the levels 1 o 3 represented by Eq.(94c) possesses 
single time of collapse and revival given by and t^_y respectively. 

Proceeding in a similar way, the population inversion for Case-II are given by 

^ 8 8 ^ V2n + 1 V2n+1 

Wg,o)it) = 7 + ^e-^{^— [7fe]> X cos[M^ +nsin[^j£^]], (976) 



W^Sfo)W = + ie-^{^--[7fcl} X cos[?4^ +nsin[^l^] 
In this case we see that the collapse and revival time arc given by 



(97c) 



V^2(2n + 1) _ 2 
5' 



2V2n + 1 2V2^ 
*U0) = ^ - (98^) 

Finally for Case-Ill, which is similar to Case-I, the population inversions among various levels are, 

H^H (,) ^ 1 _ 1 -Mi-cos[^]} ^ ,^,^9in+^ 
^^^+>^ ' 8 2 ^^/2IT3 + 

3 _2fi{i-cos[-24=]} f25(n + 3)i , _ . f 2gt 

(i) = _1 1 -Mi-cos[^]} ^ ,,,[5^n + 3)t 
^^(+)^ ^ 8 2 ^\/2^rr3 ^x/2^:T3^^ 

3 -U-co.[^]> , eosl^^l^ + .sin[^]], (996) 
8 V2n + 3 \/n — 1 

WE,,,it) = _e-M--,^]> ^ eos[^^|^ + ^sin[^=|L^]], (99c) 
^ ' V2n-|-3 v2n-|-3 

and two times of collapse are found to be 



^/2{2n + -i) _ 2 



= y^y^i;^ ^ (looa) 



_ tl^ _ VM + 3 ^ 1 



with the corresponding revival times, 

2V2n + 3 2\/2^ 



tu^\ = TT ~ TT, (101a) 

^"^^ 9 9 



i2(+) = ^ = ^ - (101&) 



respectively. Thus in the high field approximation, two distinct collapse and revival situations with 
two different time scales are present in the numerical studies also. 
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6. Conclusion 



In this paper wc have reconsidered the exact solution for ah three scmiclassical and quantized 
three-level systems where the Hamiltonians are written in the SU{3) basis. The semiclassical 
systems are solved by using the Bloch equation approach and it is shown that at zero detuning, the 
eight dimensional Bloch sphere is broken into two disjoint sectors due to the existence of two distinct 
nonlinear constants. We emphasize that such pair of nonlinear constants, which are different 
for different configurations, were not considered in earlier treatments where the Hamiltonians no 
longer obeyed a unique energy condition demanded by the SU{3) symmetry [16-20]. A significant 
outcome of our approach is that we have given a possible representation of the qutrit wavefunction 
which not only satisfies the pure state condition, but also has established its equivalence with 
the thrcc-lcvcl system. Taking the cavity modes to be in the coherent state, an Euler matrix 
based dressed state formalism is developed for finding the amplitudes of the quantized three-level 
systems and the atom-field entanglement scenario is studied by comparing the atomic entropies 
with corresponding population inversions for various initial conditions. We have explicitly derived 
the envelope functions of the equidistant cascade system for all cases and have shown that there are 
two distinct patterns of collapse and revival. The collapse and revival times of both of these patterns 
arc calculated in the high field approximation. Although the conventional quantum information 
processing has made significant progress with qubit as the primitive entity, a study of the qutrit 
wavefunction is a necessary prerequisite for developing the quantum information theory beyond 
the qubit. More exploration in this direction may give some interesting results which are within 
reach of future cavity experiments. 
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Appendix 



For completeness, in this appendix we shall present the Euler matrix based drc;ssed state for- 
malism to calculate the time-dependent amplitudes [38,39]. To find the atom-field wavefunction 
in the interaction picture, let us note that there are three possibilities for defining the initial state 
wavefunction with the quantum numbers (n, m) shown in Fig.la,b,c for the lambda, vee and cas- 
cade configuration types, respectively. Once this quantum number is assigned to any one particular 
level, then the photon quantum number of other two levels will be determined by the selection 
rules I An I = 1 and |Am| = 1, respectively. Thus at t = there are three possible choices for the 
three initial wavefunctions: (i) when the system is initially in the lower level. 



V'af(O) = 10 >(^J2J2 CmCn\m, n >, (A.la) 



771=0 71=0 

(ii) when the system is initially in the middle level. 



771=0 71=0 

(iii) when the system is initially in the upper level, 



iP+^{0) = C^\+ >^J2J2 CmCn\m,n > . (Ale) 



m— n— 

Taking above states as the initial states, the instantaneous wave function for the lambda, vee and 
cascade models can be obtained. Below we shall outline the derivation of the amplitudes of the 
atom-field wavefunction of the lambda system in detail and the remaining two are similar. 

(a) The lambda system: 
In the interaction picture, the atom-field wavefunctions for the lambda system are given by, 

OO OO C30 



i^AFit) = C- E {E CmCnd-%{t)\m,n,-> + E CmCnd-\^+,{t)\m-l,n+ 1,0 > 

n=0 m=Q m=l 

CSO 

+ E CmCnd;;,\„{t)\m-l,n,+ >}, (A2a) 

V-AfW = Co^ E i E CmCnd°Al,n-l{t)\m + 1, n - 1, > +CmCn E Wl"^' + > 

771=0 n=l n=0 

OO 

+ J2<^mCnd°J^,n-Mm,n -!,+ >}, (A26) 



71=1 



^Ipit) = C+ E {CmCnd+'l,Jt)\m + l,n,-> +CmCnd+%+,{t)\m,n + 1,0 > 

m=0,n=0 

+CmCnd+%it)\m, n, + >}, {A.2c) 

where (t) etc represent time-dependent amplitudes which are to be evaluated for three possible 
choices of basis states shown in Fig. la. Using superposition principle, the total wave function of 
the atom-field system is given by 

Krit) = rA^t) + i,%{t) + Vl^(t). (A.3) 

Finally, rearranging the summation of the eigen modes, the atom-field wavefunction is given by 



22 



E 

n,m=0 



{C-'l{t)\n,m,- > +C°^\{t)\n,m,0 > +C+:\it)\n,m, + >}, 



(AA) 



where C+^(t) etc are the time-dependent ampUtudes in the atom- field entangled basis given by 

C'm^(i) = C^CmCnd^ni^) + CQCm-lCn+ldPm-l,n+l{^) + C'+C'jn-lC'n<^m-l,n(^)' (A.So) 

To derive above expressions, wc note that the interaction Hamiltonian of the lambda system in 
the number state basis is given by {H\^)m,n =< 4'^af\^i'\''Pa'f > '^^^ be diagonalized as 



{A.6) 



where Af]^,A^ and (i = — ,0, +) be the cigcn values of the Hamiltonian and „ bo the 
transformation matrix which gives the dressed wave function from the bare states by the orthogonal 
rotation, 



(A.7) 



|m 


n,l> ' 




m — 1, n, — > 


1 rn 


n,2 > 


T^i 


|m - l,n-h 1,0 > 


_ \m 


n, 3 > 




\m,n, + > 



In Eq.(A.7), the orthogonal matrix may be parameterized by the Euler matrix 



a 



11 

5^21 



^31 



a 



12 

^2 



*32 



a 



13 

43 



^33 



(A.8) 



where. 



'^ii 

^2 



a. 



*13 
21 



a 



'^31 
i 

32 



-l'-2 
-1*2 



'-3^2^1 
■ C3C2S1 



"22 = 
"23 = 

a" 



13 

-Sl4 - 4«2Cl 

— S-1 Sn + CoCnC-i 



(A.9) 



-1^3 

*3*2 



'3"-2 



= C3 



with = sin 6'^ and = cos^^ {p = 1,2,3) three mixing angles which are to be evaluated 
for the above choices of the basis states shown in Fig.l. For choice-I, the time-dependent proba- 
bility amplitudes in Eq.(A.2a) obtained from the initial wavefunction Eq.(A.la) are given by the 
transformation. 



d-m-lAii^) 
dm-l,n+li^) 

dMt) 













. (A.lOa) 



e 

Similarly, for choice-II in Eq.(A.lb), the corresponding amplitudes in Eq.(A.2b) are obtained by 
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e-*^03t 

e"*^02t Q 



e^^oi* 

and finally for choice-Ill, the amplitudes in Eq.(A.2c) are obtained from Eq.(A.lc) 



(A.lOb) 















' 1 ' 


m,n 










T+ 





















(A.lOc) 



Finally plugging Eqs.(A.lO) back into Eq.(A.2) we obtain required time-dependent amplitudes of 
the lambda system given by Eq.(A.5). 

(b) The vee system: 
Proceeding in similar way, the atom-field wavefunction for the vee system is given by, 

oo 

^AF{t)= E {^?„-j;(t)Krn,->+C'0;;,(t)Km,0>+C'+^(t)|n,m,+ >}, (All) 

n,m=0 



where the time-dependent amphtudes (t) in the atom-field entangled basis are given by 



CZit) = CYCmCr^d-lit) + Co^C„C„_i<„_i(t) + ClCm-iCr.d+l.Jt), {A.12a) 

^^,nit) = C- CmCn+ldjy^n+li^) + (^0 ^mCnd'^nit) + C+Cm-lCn+ld^^ ^_^_i{t), {A.12b) 

C^m^nii) = CY C'm+lCnd^i,^{t) + Cq Cm+lCn-ld^m+l,n-l{t) + C]. C'mCnd^n{t) . (>1.12c) 

(c) The Cciscade system: 

Similarly, for the equidistant cascade system the wavefunction is found to be, 

oo 

^AF{t) = E{C--(f)|n,- > +Cf (t)|n,0 > +C+^{t)\n,+ >}, (A13) 

where the time-dependent amplitudes are given by 

C-^{t) = ClCnd-~{t) + Cfa+id?,^i(t) + ClC^+2d+Ut), {A.lAa) 

Cf (t) = C=C„_irf-=i(t) + Co^C„df (t) + ClCn+idllS), (A.146) 

Ct^{t) = CtCn-2d-%{t) + Cf C„_irfO= + ClCnd-irit). {A.1AC) 

Now the primary task is to calculate the values of (i^„(t) for the different Hamiltonians of the 
three-level systems and hence to obtain the corresponding time-dependent amplitudes. 
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FIGURE CAPTION 



Fig.l: The energy levels of the lambda, vee and cascade transition types with the 
energy levels arranged as E3 > E2 > Ei. For each configuration, three choices of the 
basis states are possible (shown for lambda system only.) 

Fig. 2: The time evolution of the atomic entropy and population inversion of the 
lambda system for Ccise-I (Red), Case-II (Green) and Case-Ill (Blue) with 

gi = -2, 92 = -1, Oim = \/30 and a„ = \/26. 

Fig. 3: The time evolution of the atomic entropy and population inversion of the vee 
system for Case-I (Red), Case-II (Green) and Case-Ill (Blue) with the same values 

of 91, 92, ctm and q;„ as for Fig. 2. 

Fig. 4: The time evolution of the atomic entropy and population inversion of the 
equidistant cascade system for Case-I with g = .1 and Q!„ = %/35- The presence of two 
distinct type of collapse and revival cu-e evident for W12 and Wfa transitions only. 

Fig.5: The time evolution of the atomic entropy and population inversion of the 
equidistant cascade system for Case-II for the same values of g and a„ as for Fig. 4. 
Here we have single collapse and a single revival time. 

Fig. 6: The time evolution of the atomic entropy and population inversion of the 
equidistant cascade system for Case-Ill with the same values of g and Q!„ as for 
Fig. 4. The plots are similar to those for Case-I. 
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^ h{uJi + W2) 



E^ ^ -hioi 




Choice-I Choice-II Choice-Ill 

|m— + > |m, n — !,+ > \n,m,+ > 



|m— l,,n + l,0> \m,n,0> 



|m, n + 1, > 



|m, n, — > 



|m + 1, n — 1, — > |m + l,n,— > 



Fig. la 



E^ — huJi 



EY = hjJr> 



EY = -(wi +W2) 













Til 



Ef — hu!2 

E^ = h{uji - UJ2) 
Ef = —huji 



Fig. lb, c 
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